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Abstract 

This paper is concerned with the derivation and analysis of hydrodynamic models 
for systems of self-propelled particles subject to alignment interaction and attraction- 
repulsion. The starting point is the kinetic model considered in [10] with the addi- 
tion of an attraction-repulsion interaction potential. Introducing different scalings 
than in [10], the non-local effects of the alignment and attraction-repulsion interac- 
tions can be kept in the hydrodynamic limit and result in extra pressure, viscosity 
terms and capillary force. The systems are shown to be symmetrizable hyperbolic 
systems with viscosity terms. A local-in-time existence result is proved in the 2D 
case for the viscous model and in the 3D case for the inviscid model. The proof 
relies on the energy method. 
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1 Introduction 

The context of this paper is the hydrodynamic limit of a kinetic model for self-propelled 
particles. The self-propulsion speed is supposed to be constant and identical for all the 
particles. Therefore, the velocity variable reduces to its orientation. The particle interac- 
tions consist in two parts: an alignment rule which tends to relax the particle velocity to 
the local average orientation and an attraction-repulsion rule which makes the particles 
move closer or farther away from each other. This model is inspired both by the Vicsek 
model [28] and the Couzin model [2, 8]. 

The model studied in this paper is a generalization of the model of [10] with the addi- 
tion of an attraction-repulsion interaction potential. More importantly, a different scaling 
is investigated. In this scaling, the non-local effects of the alignment and attraction- 
repulsion interactions are kept in the hydrodynamic limit and result in extra pressure and 
viscosity terms. Beyond the statement of the model, the main result of the present paper 
is a local-in-time existence theorem in the 2D case for the viscous model (when the non- 
local effects are retained) and in the 3D case for the inviscid model (when the non-local 
effects are omitted). Both proofs rely on a suitable symmetrization of the system and on 
the energy method. 

There has been an intense literature about the modeling of interactions between in- 
dividuals among animal societies such as fish schools, bird flocks, herds of mammalians, 
etc. We refer e.g. to [1, 2, 8, 17] but an exhaustive bibliography is out of reach. Among 
these models, the Vicsek model [28] has received particular attention due to its sim- 
plicity and the universality of its qualitative features. This model is a discrete particle 
model (or 'Individual-Based Model' or 'Agent-Based model') which consists of a time- 
discretized set of Ordinary Differential Equations for the particle positions and velocities. 
A time-continuous version of this model and its kinetic formulation are available in [10]. 
A rigorous derivation of this kinetic model from the time-continuous Vicsek model can be 
found in [3] . In the present paper, we extend this model by adding an attraction-repulsion 
force. 

Hydrodynamic models are attractive over particle ones due to their computational 
efficiency For this reason, many such models have been proposed in the literature [5, 6, 
7, 14, 19, 20, 26, 27]. However, most of them are phenomenological. [10] proposes one 
of the first rigorous derivations of a hydrodynamic version of the Vicsek model (see also 
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[18, 23, 24] for phenomenological derivations). It has been expanded in [11] to account for 
a model of fish behavior where particles interact through curvature control, and in [12] to 
include diffusive corrections. Other variants have also been investigated. For instance, [15] 
studies the influence of a vision angle and of the dependency of the alignment frequency 
upon the local density. [9, 16] propose a modification of the model which results in phase 
transitions from disordered to ordered equilibria as the density increases and reaches a 
threshold, in a way similar to polymer models [13, 21]. 

The organization of the paper is as follows. In section 2, we introduce the model of 
self-propelled particles and set up the associated kinetic equation. We then discuss various 
scalings which lead to the derivation of the studied hydrodynamic models. We introduce 
four dimensionless parameters in the problem: the scaled interaction mean-free path s, 
the radius of the interaction region t], the noise intensity 5/e and the relative strength 
between the attraction-repulsion and the alignment forces s. The scaling considered in 
[10] ignores the attraction-repulsion force and supposes that e = rj — > 0, 8 = 0(1). Here, 
we investigate four different scaling relations. 

1. The weakly non-local interaction scaling without noise: r\ = y/e, 5 = 0, s = rj 2 . 
The resulting model is a viscous hydrodynamic model with constrained velocity on 
the unit sphere. For this case, we assume that the solutions of the kinetic equation 
are monokinetic. We justify this assumption by studying the space homogeneous 
kinetic model and prove that the solutions converge on the fast e time scale to the 
monokinetic distribution. We also highlight the variational structure of this space 
homogeneous kinetic model. Note that the scaling assumption rj = yfe is different 
from the one used in [10]. It corresponds to increasing the size of the interaction 
region in the microscopic variables by a factor 1/y/e, as e — > 0. Therefore, more and 
more non-local effects are picked up in the hydrodynamic limit. These non-local 
effects give rise to the viscosity terms in the macroscopic models which make an 
original addition from previous work. 

2. The local interaction scaling with noise. This is the scaling proposed in [10] which 
is recalled here just for the sake of comparisons. It consists in r, <C e, 5 = 0(1), 
s < rj 2 . The resulting model is the inviscid hydrodynamic model with constrained 
velocity on the unit sphere. 

3. The weakly non-local interaction scaling with noise. This scaling unifies the two 
previous scalings. It consists in rj = yfe, 5 = 0(1), s = rj 2 . Again, the resulting 
model is a viscous hydrodynamic model with constrained velocity on the unit sphere, 
but with modified coefficients as compared to the first scaling. We note however, 
that in the zero noise limit 5 — > 0, we recover the system obtained with the first 
scaling, which provides another justification of the monokinetic assumption in the 
derivation of the model. 

4. Capillary force scaling. This corresponds to rj — y/e, 5 = 0(1), s — 1. Therefore, 
here, the attraction repulsion force is of the same order as the alignment force. 
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However, we make the additional assumption that the zero-th order moment of the 
potential is zero, which expresses some kind of balance between the attraction and 
repulsion effects. This results in a model like in the previous scaling, but with the 
addition of a term analog to the capillary force, induced from the attractive part of 
the potential. 

In section 3, we prove local well-posedness for all the models derived in section 2, 
except the last one (capillary force scaling). All the remaining systems have the same 
form of a symmetrizable hyperbolic system with additional viscosity. In section 3.1, we 
prove the local-in-time existence of solutions for the viscous system in 2D and in section 
3.2, we show the same result for the inviscid system in 3D based on the energy method. 
Finally, a conclusion is drawn in section 4. 



2 Derivation of hydrodynamic models 

2.1 Individual-Based Model of self-alignment with attraction- 
repulsion 

The starting point of this study is an Individual-Based Model of particles interacting 
through self-alignment [28] and attraction-repulsion [2, 8]. Specifically, we consider N 
particles Xk € M d moving at a constant speed v k £ S '" 1 . Each particle adjusts its velocity 
to align with its neighbors and to get closer or further away. The evolution of each particle 
is modeled by the following dynamics: 

— = v k (2.1) 
dt y ! 

dv k = P v ±(y k dt + V2ddB*). (2.2) 
Here, P v ± is the projection matrix onto the normal plane to v k : 

P v ± = Id — v <g) v. 

It ensures that Vk stays of norm 1. B\ is a Brownian motion and d represents the noise 
intensity. Both the alignment and attraction- repulsion rules are encoded in the vector v k : 

jk + r k 

Vk 



\jk + r k y 

where j k counts for the alignment and r k for the attraction-repulsion: 

j k = ^K[\xj-x k \)vj , r k = Ys®{\xj-x k \) * i ~_ Xk (2.3) 

\Xj Xk\ 

J 3 

The kernel K is a positive function, $' can be both negative (repulsion) and positive 
(attraction). In figure 1, we give an example of functions K and modeling the popular 
"zone-based" model for fish behavior [2, 8, 22]. 
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K(r) 
&( r ) 



repulsion ■ alignment 



attraction 



Figure 1: An example of functions K and $' modeling the "zone-based" model: repulsion- 
alignment-attraction. 



When the number of particles becomes large (i.e. N —> oo), one can formally derive 
the equation satisfied by the particle distribution function f(x,v,t) (i.e. the probability 
distribution of the particles in phase-space (x,v)). Under suitable assumptions [3, 10, 25], 
/ satisfies: 



ft + v-VJ= -V v ■ [{P v xv f )f] + dA v f, (2.4) 



where 



Jf + r f 



\jf + r f\ 

j f = I K(\x' - x\)v' f(x',v',t)dx'dv' 



r f = -Va, / $(\x' -x\)f(x',v',t)dx'dv'. 

Jx',v' 



The function $ is the antiderivative of $' which vanishes at infinity (i.e. $(r) 0). 
Using the distribution /, we want to identify the asymptotic behavior of the model in 
different regimes. This is the purpose of the next section. 

2.2 Scaling parameters 

Introducing two dimensionless parameters e and t], the starting point is the following 
scaled version of the previous kinetic model for the distribution function f(x, v, t): 

ft + v- V x f = --V„ • [(P v ±v])f] + dA<,f, (2.5) 
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where 



v f - 






■n 

3f = 


J 

J(x>, 


r) = 





2 r V 

<nV?T' ( 2 ' 6 ' ) 
f +T) r f\ 



\x' — X\ 



K ± L ) v' f(x', v', t) dx' dv\ (2.7) 

jn x§ n-l \ T) 



1 

2n 



$ ^ ^ /(a/, v', t) dx' dv'. (2.8) 

(ac',t;')eR"xS n - 1 \ V J 

The first term (given by jl) expresses the alignment interaction (like in the Vicsek dy- 
namics [28]) while the second term (given by rj) expresses the repulsion interaction (like 
in e.g. [2]). The expression of means that the alignment term will prevail in the limit 
77 — >■ 0. In this paper, we will consider various possible assumptions concerning the relative 
speeds of convergence of e and 77 to 0. 
We denote: 

K(\£\)d{ = k , 

We can always assume that k = 1. The potential $ is said to be repulsive if <f> > 0. 
Defining the moments pf and PfUf of / by 

P f = f(v')dv', 

pju f = / f(v')v'dv', 

we have the following Taylor expansion of vj: 

v} = n f + v 2 ^-£ f + o( v % 
r pf\ u f\ 

n f = £ f = P a ±(kA{p f u f ) - 4>V xPf ). 

Inserting this expression into the kinetic equation (2.5), we get 
/ t + «-V t /=~V,-[(P^/)/] 

2 1 2 

p-rV^ • [{P^if)f\ + dA v f + o(5_). (2.9) 

e pf\u f \ e 

We now consider three different scaling limits which lead to models for which we will 
prove local existence of classical solutions. 
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2.3 Weakly non-local interaction scaling without noise 

In this scaling limit, we assume no noise d = and the following ordering between the 
two parameters e and rj: 

2 

e^O, -j -»■ 1. 

/ e satisfies (keeping only the 0(1) terms in e or larger): 

ft + u ■ V x f = --V v ■ [(P v ^ f e)r] ^V t , • [(P v x£ f e)f] . (2.10) 

For / e to converge, we need to assume that the leading order term at the right-hand 
side of (2.10) vanishes, i.e. that f e satisfies: 

V v -[(P v M fe )f} = 0, \/e>0. 

This is equivalent to assuming that f e is a monokinetic distribution, i.e. 

f £ (x,v,t) = p £ (x,t)5(v,n £ (x,t)), (2.11) 

where S(v,v) is the delta distribution on the sphere at the point v. The assumption of 
monokinetic distribution requires, in order to be consistent, that there is no noise, which 
is the reason for assuming d = 0. 

Proposition 2.1 For monokinetic solutions (2.11), p and Q are independent of e and 
satisfy the following system: 

d tP + V x -( P n) = 0, (2.12) 
dt(pn) + V, • (pfi ® fi) + (j)P^VxP = kP u xA(pu). (2.13) 



Proof. The result follows from multiplying (2.10) by 1 and v and using Green's formula. 



Remark 2.1 The repulsive force contributes for a pressure term at the left-hand side of 
the momentum equation, which otherwise would not be strictly hyperbolic, and would fall 
in the class of Pressureless Gas Dynamics models [4]. 

In order to justify the monokinetic assumption, we consider the spatial homogeneous 
equation: 

ft = ~ ~V„ • [(P„xfy)/] , (2.14) 

and show that its solution relaxes to a monokinetic distribution (2.11) at the fast e time 
scale. More precisely, we have 
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Proposition 2.2 We assume that J f\t=odv = 1 so that we have j f(v,t)dv = 1 for all 
times. Therefore, u = J v f(v)dv (we omit the index f when the context is clear). We 
also assume that u t=0 ^ 0, otherwise, the dynamics is not defined. Then, f{t) — > /(oo) 
where /(oo) is of the form /(oo) = S(v, Q). 

Proof. We introduce the variance: 

Hf) = I \v-u\ 2 f(v)dv 



\u 


2 _ 


(v ■ u) 2 




u 





(u ■ v)) f{v) dv. 

We have J-(f) = 1 — \u\ 2 . We note that \u\ 2 is the classical order parameter [28]. In 
Lemma 2.3 below, we prove that Hf) satisfies the following dissipation equation: 

Hf)t + - / H2 ~ ( | ' U? f(v) dv = 0. (2.15) 
Since \u\ 2 — (v ■ u) 2 > 0, we have 

Hf)t < o. 

Therefore, F{f) is a decreasing function of time. Furthermore, if / is a distribution such 
that 

f(v)dv = 0, 

then, either / = 5(v, Q) or / is of the form: 

/ = a5(v,n) + (l-a)5(v,-Q), 

where a G [0, 1] and Q G S n_1 . The second form is called a dipole. So, unless / is a delta 
or a dipole, F{f) is strictly decaying. Since a dipole is unstable, it is never reached in the 
course of the dynamics. Therefore, Hf) — ^ as t — )• oo. Since Hf) = is equivalent to 
/ = 5(v, Q), this shows that / — > 5(v, Q) as t — > oo and the typical convergence time is s. 



Lemma 2.3 Any solution f of (2.14) satisfies (2.15). 
Proof. We write 

sHf)t = f(l-(u- v)) d t f{v) dv-Uf J v f(v) dv 

(1 — (u ■ v)) d t f{v) dv — J v d t f(v) dv ■ u 
d t (J f(v) dv)-2 J(u-v)d t f(v) dv 
-2 J (u-v)d t f(v)dv. 
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In the last equality, we have used that f f(v,t)dv = 1. Now, multiplying (2.14) by 
— 2(u ■ v), integrating with respect to v and using Green's formula, we get: 



But 



u-P v ±Q = rr(M 2 - (u-v) 2 ), 



which leads to the result. 



is a free energy for the problem (2.14) and provides a variational structure. First, 
let us denote by v = V/-J 7 the gradient of T with respect to /. It is defined by 

where g is an increment of /, i.e. a function g(v) satisfying J g(v) dv — (so that / + g 
satisfies the admissibility condition f (/ + g) dv — 1). 

Proposition 2.4 Eq. (2.14) can be recast as 

ft - -^V, • {\V v (V f F)]f} = 0, (2.16) 

which shows that the flow of (2.14) has a gradient flow structure in the Wasserstein 
metrics [30]. We have: 

Hf)*+^7 / |V,(V / ^)| 2 /^ = 0. (2.17) 



Remark 2.2 Eq. (2.17) provides another proof of the decay o/J r (/) with time. 
Proof. We compute: 

TT^^ -2 ^ / g(v)vdv = (-2(u ■ v),g), 



Sf 

which yields 

v = V/J 7 = —2u ■ v. 

A simple computation shows that 

\7 v v = —P v ±u = — \u\P v ±Q. 

Therefore, eq. (2.14) can be written as (2.16). Now, multiplying by V/J 7 , integrating 
over v and using Green's formula, we get (2.17), which ends the proof. i 
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2.4 Local interaction scaling with noise 

In this scaling we assume that d = 0{j). More precisely, we let: 

d — -, 
e 

with 5 a given constant. We also assume that e and r\ are such that: 

2 

e^O, n-)-0, — ^0. 

e 



With this last assumption, the 0(— ) term in (2.9), which results from the non-locality 
of the average alignment direction, vanishes. Therefore, this scaling keeps only the local 
contribution of the alignment interaction. The resulting asymptotic problem, keeping 
only terms of order 0(1) or larger, is written: 

ft + v ■ VJ e = ~\ {V, ■ [(iVa)f ] + 5A v f} . (2.18) 

The limit of (2.18) as e — > has been studied in [10] in dimension 3 and in [15] in any 
dimensions. The result is stated in the following theorem. 

Theorem 2.5 We have f e — > pM^ where Mq(v) is the Von Mises-Fischer distribution: 

, i exr>(8(v ■ Q)) dv 1 . 

M n (v) = - )J , /3 = -, (2.19) 

J v& n-i exp(/3(t; • n)) dv 5 

and p and Q satisfy the following system: 

dtp + ci V x ■ (pQ) = 0, (2.20) 

P (d t n + c 2 vt ■ v x n) + sp a xV xP = o. (2.21) 

The constants c\ and c 2 are defined by 



c y = / M n (v)(v-n)dv, 

f v&n ^ M n (v) h(v ■ O) (1 - (v ■ n) 2 ) (v ■ Q) dv 



c 2 



M a {v) h(v ■ O) (1 - (v nf) dv 



where h(v-Q) is the Generalized Collision Invariants (GCI) [10] and is defined as follows 
in the n- dimensional case [15]. Set ip a (v) = h(Q ■ v ) (a • v ) where a G M n is any vector 
such that a ■ Q = 0. Then, ip a is the unique solution in the Sobolev space if 1 (S n ~ 1 ) with 
zero mean, of the following elliptic problem: 

-A v i/) - (3(Q ■ V„)V> = a ■ v . 



Proof. We refer to [10] in the three dimensional case and [15] in the general n-dimensional 
case. ■ 
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2.5 Weakly non-local interaction scaling with noise 

In this section, we propose a scaling which unifies the two previous ones. In this scaling 
we assume that d — | , with 5 is a given 0(1) constant and that: 

2 

£^0, 77->-0, y ->■ 1. 

2 

Here ^ — )■ 1 instead of like in the previous section. Inserting these assumptions into 
(2.9), and keeping terms of order 0(1) or larger, we get 

ft + v- V x f + -^^v„ ■ [(P^/O/l = 

Pfe\Ufe\ 

- 1 - {V, ■ [(/VfyOf ] + SA v f} . (2.22) 
The limit e — > is investigated in the following theorem: 

Theorem 2.6 We have f e —> pM^ where Mq(v) is the Von Mises-Fischer distribution 
(2.19). p and Vt satisfy the following system: 

dtp + cxV x ■ (pO) = 0, (2.23) 

P (d t n + c 2 n ■ v x n) + (5 + c^)p^w xP = c 3 k Cl p^A( P n), (2.24) 

where the constants c\ and c 2 are defined as in Theorem 2.5 and 

c 3 = (n ~ 1)(5 + C2 . (2.25) 

Cl 



Remark 2.3 We notice that Ck — > 1 as 5 — > for k = 1, 2, 3 and we recover the noiseless 
system (2.12), (2.13) when 5^-0. 

Proof. We write (2.22) as 

{T l + T 2 )r = -QU £ ), 

where T\ + T 2 and Q are respectively the operators appearing at the left and right hand 
sides of (2.22). T\ = d t + v ■ V x and T 2 is the remaining part of the left-hand side. 
Integrating over v and letting f £ — > pM^ leads to the mass conservation equation (2.20) 
unchanged, since T 2 is in divergence form and vanishes through integration with respect 
to v . 

Now, to get the momentum equation, we proceed like in [10]. From the Generalized 
Collision Invariant property [10], it follows that 

P n ± / T(pM n )hvdv = 0. 
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Now, the term 

P l :=P n x / TtipM^hvdv, 



gives rise to the same expression as in [10]. This expression is 

Pi = Papd t Q + 7 fi • V X Q + aP n ±V x p, 

with 

a = — ^- / M a {v) h(l-(v n) 2 ) dv, 

1= ( * / M n (v)h(l-(vQ) 2 )(vQ)dv. 
[n - 1)6 J veS n-i 

Dividing by a/3, we find the coefficients c 2 and 5 of (2.21) (we recall that 05 = 1). 
We introduce the notation 

I ■= tpMn = iV(/cCiA(pfi) - 4>V xP ), 

and consider 



P 2 := P^ [ T 2 (pM n )hvdv 

-P^ [ V v -[(P v x£)M n ] hvdv. 



Using Green's formula, we get 



P 2 = --iV / [P^-V^/i^Mndt;. 
We note that (P v ±£) ■ V v <fr = £ • V„</>, with being any component of hv. We deduce that 



P 2 = / (£-V v )(hv)M n dv 

T 



Cl 

Now, we use the formulas: 



-P^ f / V w (/it;)Mndi; > ) 1 (2.26) 



/ V v gdv — (n — 1) / i^cfo 

/ (V v g)hdv — (n — 1) / vghdv— / (V v h)gdv 
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for any pair of scalar functions g, h on S n 1 . We recall that V„Mn = PP v ±QMq. Since 
£ ■ Q = 0, we compute the matrix 



D := ( f V v (hv)M a dv) P a . 

= (n-l)(/ («®v)/iM n du ) Pqj 



-/3 ( / (v ® P w x.fi)/i M n tfo J P n x 

■- ( n - 1)D X - D 2 . 
We decompose 

f = + f j_ = P^xf , U|| = (u • 

Using this decomposition and the fact that integrals of odd degree polynomials of v± over 
S n_1 vanish, we have: 

£>!=(/ ( u± ® v± )hM Q dv) P n ± = aP n ±, 

and 

D 2 = p( [ {(v\\+v x )®(£l-{n-v)){v\\+v ± ))hM sl dv)p fl x. 

Owing to the fact that any term of the form (A <g) v\\)Pq± = for any vector A, we have 
since v\\ is parallel to Q: 



D 2 = (v± ® v ± )(Q ■ v) hMndv) P^x = - 7 P nJ 

Inserting these results into (2.26), we get 

(n — l)a + 7 



P 



2 



Cl 



Collecting all the results and dividing by a/3, we are led to the momentum equation (2.24), 
which ends the proof. i 



2.6 Attraction-repulsion potential: induced capillary force 

In this section, we investigate the case where the attraction-repulsion force term is of the 
same order as the alignment term in the expression of the alignment direction v 7 }, i.e. we 
assume that 

\jf + r f\ 
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where jj and rj are respectively given by (2.7) and (2.8). Note that, by contrast to (2.6), 
there is no rf in front of in (2.27). 

The Taylor expansion of v^j is now given by 

fl/ ^r!^, , * / = P^(*A(p / u / )-&V a ,Ap / ), 



where 



h 1 



.r' 



Here, we suppose like in [2], that the potential is repulsive at short scales and attractive 
at large scales (See Fig. 2). Therefore, $(|£|) is supposed to decrease for |£| G [0, £*] and 
to increase for |£| G [£*, +oo). Furthermore, since $(|£|) is supposed integrable on IR n , we 
have $(|£|) ^ as |f| oo. It results that < and that $ > for |f| G [0,f ] 

and $ < for |£| G [£o, +oo) where £ < £*• We make the additional assumption that the 
zero-th order moment vanishes: 

= 0, 

which expresses the balance between the attractive and repulsive parts of $. Given the 
above assumptions, the second moment is negative: 

02 < 0. 
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With these assumptions, the Taylor expansion of vj simplifies and becomes: 

l,? = !W ^' + 0(,2) ' 

n ' = j^' (f = Pnf(kA(p,u,)-<l>2V I Ap f ). (2.28) 
Now, we can develop the same theory as before, assuming that 

2 

e^O, r/->0, y -> 1. 

Inserting these assumptions into (2.9), and keeping terms of order 0(1) or larger, we get 
(2.22) but with t given by (2.28). The limit e — > can be performed like in section 2.5 
and we obtain the following theorem: 

Theorem 2.7 We have f £ — > pMn where Mq(v) is the Von Mises-Fischer distribution 
(2.19). p and Q satisfy the following system: 

dtp + c x V x • (pSl) = 0, (2.29) 

P (d t n + c 2 n ■ v x n) + 6P Q xV xP = c 3 k Cl p n ^A(pQ) + c 3 \<p 2 \p n ^ x Ap, (2.30) 

where the constants c\ and c 2 are defined as in Theorem 2. 5 and c 3 as in Theorem 2. 6. 

Remark 2.4 The last term at the right-hand side of (2.30) has the same expression as 
the capillary force in fluid dynamics, except for the projection operator P n ± . This capillary 
force is induced from the attractive part of the potential $. 

3 Existence theory 

3.1 Existence in 2D with viscosity 

This section is concerned with a local existence result in 2D for a system of the general 
form 

d t p + v, • ( P n) = o, (3.i) 

p(d t n + ctt ■ V,fi) + P^V x (p(p)) = /iiVA(pft), (3.2) 

where the constants c G M and p > are given and the pressure relation p(p) satisfies 
p'{p) > 0. All systems derived in the previous section can by recast in this form, with 
a particular choice of c, p and p(p), after time rescaling, except for the last one (section 
2.6) involving the capillary force. The system is supplemented with initial data p > 
and fl such that \fl \ = 1. We assume that the domain is the square box II 2 = [0, l] 2 
with periodic boundary conditions. 
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Theorem 3.1 We assume that the initial data belong to H m (U 2 ) with m > 2. Then, 
there exists a time T > and a unique solution (p,<p) in L°°([0, T], H m (Jl 2 )) R H 1 ([0,T], 
H m ~ l (Jl 2 )) such that p remains positive. If, in addition, p > 0, then, the solution also 
belongs to L 2 ([0,T], H m+l (Ii 2 )). 

Proof. In 2D, we can set Q = (cos (p, simp). We recall that 

d t n = n x d t ip, \7 x -n = (n ± - v x )<p, p n ± = n ± ® n ± , 

with Q 1 - = (— simp, cosip). Then, we have 

A(pfi) = Ap n + 2tt ± (V x p ■ V x <p) - 2pQ \V x <p\ 2 + pQ ± Aip, 
D, 1 ■ A(pQ) = pAip + 2 (V x p ■ V x <p). 

Therefore, system (3.1), (3.2) is written: 

(a t + n-v x )p+p(n ± -v x )<p = o, (3.3) 

(d t + cQ ■ V, V + ^(Q ± - V x )p = p (Alp + 2 Vxp - Vxif ) . (3.4) 
p V P J 

Introduce p = a(p) and A(p) such that 



a y )= Vrtti \(p) = a\p)p, h(p) = 2\np. (3.5) 
P 

Then, system (3.3), (3.4) becomes: 

(dt + n- V x )p + X(p) (^ ■ V x )<p = 0, (3.6) 
(d t + ctt ■ V x )<p + A(p) (Q L ■ V x )p = p (Aip + V x h(p) ■ V x <p) . (3.7) 

From (3.6), we have the following a priori estimate (maximum principle): 



Pmin 

where 



exp(- / \\V x <p(-, s)\\ L oo(p2) ds) < p < p max exp( / || V x ip(-, s)||l°°(ii2) ds), (3.8) 
Jo Jo 



Pmin = min p (a;) , p max = max p (x) . 



We remind the following lemmas [29]: 
Lemma 3.2 For any pair of functions f , g in H m (W l ) fl L°°(IR n ) ; we have: 

\\fg\\H- < C (\\f\\ Hm \\g\\ L oo + \\f\\ L °°\\g\\H™) . 

If additionally, we suppose that V/ G L°°(IR n ) ; we have, for any a G N" ; with \a\ 

En 
i=i on = m : 

\\D a (fg) - fD a g\\ Hm < C(||/|WMU« + \\V x f\\ L ~ \\g\\ H ^) , 



where D a = d x <*i x ° 
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Now, with \a\ < m, we take the D a derivative of (3.6) and multiply it by D a p and 
integrate it with respect to x. Similarly, we take the D a derivative of (3.7) and multiply 
it by D a ip and integrate it with respect to x. We sum up the resulting identities. Using 
the notation 

(f,g) = / fgdx, 

we find: 



= (D a p,D a p t ) + (D a (p, D a (p t ) 

+ (D a p, D a ((n ■ V x )p)) + c{D a cp, D a ((n ■ V x )cp)} 

+{D a p, D a (\(p) (n L ■ v x )<p)) + D a (\(p) (n 1 - ■ v,)p)> 

—p,(D a if, D a Aip) 
-p>(D a <p : D a (V x h(p)-V x ip)) 
= h + ... + h 

Then: 



where || • || just indicates an L 2 norm. Now, for the remaining terms, we have the following 
lemma 



Lemma 3.3 We have: 

\h\ < CfllpHwi,- + II^H^i.oo) {\\p\\ 2 Hm + H^H^m), fc = 2, 3, 
|/ 5 | < ^l|V^V|| 2 + C , (||p||^ 1 , 00 + ||^|| 2 w i,o )(||p||^ m + ||^||^), 

where C denote generic constants depending on the parameters of the problem. 

The proof of the lemma is postponed at the end. 

Adding all these terms together for all possible indices a such that |a| < m, we have, 

+C r (||p||^i, 00 + imi^loc + 1) (||p||^ + \\<f\\ 2 Hm ). 

For m > | + 1, we have 

IWI.OO + ||<^|| W l l0 o < C(||p||i/m + ||<p||ffm), 



17 



and get 



2^M| m + ll* m ) + fl|V* m < C(||p||^ + ||^||^ + l) 2 . 

Gronwall's inequality leads to the local existence of a solution (p, ip) in L°°([0, T], H m (U 2 )) 
which, if p > 0, also belongs to L 2 ([0, T], H m+1 (U 2 )) and which satisfies the a priori bound 
(3.8). To get time regularity, we directly use eqs. (3.6), (3.7), take the if m_1 norm, apply 
Lemma 3.2, and find 

ll/Mli?-™- 1 + Hi?™- 1 < C||v 9 llH m + i + CGIpllw 1 . 00 + IM^ 1 ^) (||p||#™ + 1 1 1 1 ) - 

Using the previous estimates, we deduce that (p, <p) also belongs to -£^([0, T], if m_1 (n 2 )). 
The estimates on p immediately transfer to p since a(p) is smooth and invertible for p > 0. 



Proof of Lemma 3.3. Estimate of I 5 : Using Green's formula and Cauchy-Schwartz 
inequality, we have: 

\h\ < n\\VD a <p\\\\V x h(p)-V x tp\\ Hm -i 

< C\\VD a v\\ (\\p\\h*\\V<p\\l« + \\Vp\\l°°\\<p\\h>*) 

< |||vdvh 2 + c (||p||^||v ¥ »|Ucc + IIVpiUooii^n^) 2 

< Iiivdvii 2 + c(\\p\\ 2 whoo + Mlv^) (\\p\\ 2 H ™ + 11*™)- 

The second inequality uses Lemma 3.2 and the third one uses Young's inequality. 



Estimate of I^: We write 

h = (D a p, X(p) (Q 1 ■ V x )D a v ) + (D a <p, (X(p) (fi x • V x )D a p)) 
+ (D a p, {D a (X(p) (fl 1 - ■ V» - X(p) (Q 1 - ■ V x )D a y)) 

+(D a <p, {D a (x(p) (n x ■ v x )p) - x(p) {n 1 ■ w x )D a p)) 
= J1 + J2 + J3. 

Using Green's formula, we find 

|Jx| = \(V-(X(p)n ± )D a p,D°'ip)\ 

< c(\\p\\ w i,™ + ii^iiwi.oo) (iipii^m + y\\ 2 H m)- 

Now, using Cauchy-Schwartz inequality and applying Lemma 3.2, we find that J2 and J 3 
satisfy the same inequality. 

Estimate of I 2 : The proof is similar as for I 3 and is omitted. ■ 
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3.2 Existence in 3D without viscosity 

In this section, we investigate the local existence for the inviscid problem in 3 dimensions: 

d t p + V x -( P n) = 0, (3.9) 

P (d t n + eft ■ v x ft) + p n xv x (p(p)) = o, (3.io) 

where the parameters and data have the same meaning as in section 3.1. We consider the 
system in the domain II 3 = [0, l] 3 with periodic boundary conditions. 

For this purpose, we use the spherical coordinates associated to a fixed Cartesian basis. 
In this basis, denoting by 9 e [0, tt] the latitude and <p G [0,27r] the longitude, we have 

ft = (sin 9 cos <p, sin 9 sin <p, cos 9) T , 

and we let ft# and ft^, be the derivatives of ft with respect to 9 and if. We note that 



|ft#| = 1, jft^l = sin^. 



We will use the formulas 



v x • ft = n 6 ■ v x 9 + n v ■ v x ip, 

sin 9 

(ft • v x )n = ((ft ■ v x )9) n e + ((ft ■ v,)^) %, 

ft t = VL e 9 t + ft^<pt, 

where a is an arbitrary vector. 

Introduce p and A(p) as in (3.5). Then, system (3.9), (3.10) becomes: 

p t + n-v x p + \(p)v x ■ ft = 0, 

ft* + c (ft • V ffi )ft + A(p) P,iV,/5 = 0, 



or, 



Introducing 



p t + ft • V x p + A(/5)(ft • VJ + ft v • V x <p) = 0, (3.11) 
t + c (ft • + A(p) ft e • V x p = 0, (3.12) 
sin 2 6(p t + c sin 2 9 (ft • V x )<p + A(p) ft^ • V x p = 0. (3. 13) 



U 




this system is written 

A (U)U t + Ai(C0Cr a + A 2 (U)U y + A 3 (C/)^ = 0, 
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in Cartesian coordinates x = (x,y,z), where Ak(U), k = 0, ...,4 are all symmetric 
matrices and 



If sin 6* > 0, then this system is a symmetrizable hyperbolic system. We can apply 
proposition 2.1 p. 425 of [29] and the following theorem follows immediately: 

Theorem 3.4 We assume that the initial data (po,0o,(po) belong to if m (II 3 ) with m > 
5/2 with po > 0, sin 6*o > 0. Then, there exists a time T > and a unique solution 
(p,9,ip) in L°°([0,r],if m (n 3 )) n/f 1 ([0,T],/f m - 1 (n 3 )) such that p remains positive. 

4 Conclusion 

In this paper, we have derived hydrodynamic systems from kinetic models of self-propelled 
particles with alignment interaction and attraction-repulsion force. We have particularly 
focused on the inclusion of diffusion terms under the assumption of weakly non-local 
interactions. Then, we have proved the local-in-time existence of solutions for the viscous 
system in 2D and a similar result for the inviscid system in 3D. The methods rely on a 
suitable symmetrization and on the energy method. Future works in this direction will 
consist in continuing the exploration of the mathematical structure of the system and 
particularly, trying to prove local existence of the viscous system in 3D and the treatment 
of the geometric singularity near sin 6 = 0. Another direction of work will consist of 
the numerical quantification of the viscosity as a consequence of the non-locality of the 
interaction. 
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